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Abstract 

We construct new elliptic solutions of the restricted Toda chain. These solutions 
give rise to a new explicit class of orthogonal polynomials which can be considered as 
a generalization of the Stieltjes-Carlitz elliptic polynomials. The recurrence coefficients 
and the weight function of these polynomials are expressed explicitly. In the degenerated 
cases of the elliptic functions the modified Meixner polynomials and the Krall-Laguerre 
polynomials appear. 
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1. Restricted Toda chain and orthogonal polynomials 



The main purpose of the present paper is explicit construction of orthogonal polynomials which 
are a generalization of the famous Stieltjes-Carlitz orthogonal polynomials connected with el- 
liptic functions. For theory of these polynomials and history of their discovery see, e.g. [15], 
[7], [25]. 

Our main tool will be connection between orthogonal polynomials depending on an addi- 
tional "time" parameter t and solutions of the so-called restricted Toda chain. This connection 
allows us first to construct some explicit "elliptic" solutions of the Toda chain and then to 
reconstruct corresponding orthogonal polynomials and their orthogonality measure. As we will 
see, the Stieltjes-Carlitz polynomials appear to be a very special case of these constructed 
orthogonal polynomials corresponding to a "zero time" t — case. We derive also explicit re- 
currence coefficients u n (t), b n (t) for the obtained orthogonal polynomials. These coefficients are 
expressed in terms of the Weierstrass elliptic functions in n and t . Recall that the recurrence 
coefficients of the Stieltjes-Carlitz polynomials are expressed in terms of linear and quadratic 
polynomials in n. We show that obtained polynomials are orthogonal on the whole real axis 
with a positive discrete measure. The measure is constructed explicitly. It appears that the 
parameter t can take values only inside of some interval (so-called admissible interval) in order 
for the measure and the recurrence coefficients will be well defined. The Stieltjes-Carlitz case 
corresponds to the middlepoint t = of this interval. 

We recall basic definitions and results concerning relations between Toda chain and orthog- 
onal polynomials [21], [4], [17]. 

The Toda chain equations are [22] 

u n = u n (b n - b n -i) , b n = u n+1 -u n (1.1) 

with additional condition 

uo = (1.2) 

where the dot indicates the differentiation with respect to t. In what follows we will call 
equations (1.1) with restriction (1.2) the restricted Toda chain (TC) equations. 

Let P n (x;t) be orthogonal polynomials satisfying the three-term recurrence relation 

P n +i(x) + b n P n (x) + w n P„_i(a;) = xP n (x) (1.3) 

with initial conditions 

P = l, P 1 (x)=x-b . (1.4) 

We will assume that u n ^ 0, n — 1, 2, . . .. By the Favard theorem [8], there exists a nonde- 
generate linear functional a such that the polynomials P n (x) are orthogonal with respect to 
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it: 



a(P n (x)P m (x)) = h n S nm , 



[1.5) 



where h n are normalization constants. The linear functional a can be defined through its 
moments 

c n = a(x n ), n = 0,l,.... (1.6) 

It is usually assumed that cq — 1 (standard normalization condition), but we will not assume 
this condition in the followings. So we will assume that c is an arbitrary nonzero parameter. 
Introduce the Hankel determinants 



D n = det(cj +j )ij =0 ,..., n -i, D = 1, D x = c . 
Then the polynomials P n (x) can be uniquely represented as [8] 



Pn(x) = 



C Ci 

Cl c 2 

Cn— 1 C n 

1 X 



c n+l 
C2ra-1 



The normalization constants are expressed as 



y n+l 



, /i = Di = cq. 



The recurrence coefficients «„ satisfy the relation 



u 



n 7 



Thus we have 



;i.7) 



(1.9) 



;i.io) 



(1.11) 



Assume now that the polynomials P n (x; t) depend on a real parameter t through their recur- 
rence coefficients u n (t), b n (t). Then the restricted Toda chain equations (RTE) are equivalent 
to the condition 



Pn(x;t) = -u n P n ^(x;t). (1.12) 

It is possible to choose initial moment Co(t) (normalization) such that the RTE are equivalent 
to the very simple condition 

c n = c n+ i, (1-13) 
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Hence, for the Toda chain case, the Hankel determinants D n = D n (t) have the form 

D n (t) = det(4 l+fc) (t)) iifc=0i ..., n _ 1 , D = 1, D 1 = c , (1.15) 

where means the j-th derivative of co(t) with respect to t. 

Under this condition, the RTE are equivalent also to the equations 



d 2 \ogD n D n ^D n 



' ra=l,2, .... (1.16) 



dV Dl 

Note also that for the Hankel determinants of the form (1.15) we have the useful relation 

, Dn+l D n 

bn= n ~n (L17) 

or, equivalently, 

b n = K/h n . (1.18) 

In particular, for n = we have from (1.18) 

b = ^. (1.19) 

c 

The relation (1.19) allows us to restore co(t) if the recurrence coefficient bo = bo(t) is known 
explicitly from Toda chain solutions (1.1). 

The Stieltjes function F(z) is defined as a generating function of the moments [8] 

If moments c n depend on t according to the Toda Ansatz (1.13), we then have 

F(z;t) = C ^ + ^ + ... + ^ + ... = zF(z)-co. (1.21) 

z z z z n 

In fact, the relation (1.21) is equivalent to the restricted TC equations (1.13). 
We consider also a so-called ^-generating function of another type: 

$ (P) = X>^- (1-22) 
k=o K - 

The relationship between functions F(z) and $(p) is given by the (formal) Laplace transform: 

F(z) = Ckz-"- 1 = £ c k P -^dp = / e-* z Q(p)dp. (1.23) 
fc=0 k=0 Jo kl Jo 
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For the case of the RTE with condition (1.13) we see that generating function $(p) is given 
automatically by the formal Taylor expansion 

oo k oo k 

*(p; t) = £ c k {t) p - = £ $\t) p - = c (t + P ) (1.24) 

fc=0 K - k=0 K - 

of co(t +p). Thus the ^-generating function is given just by the shifted co(t +p) zero-moment 
function. The Stieltjes function is given then as the Laplace transform 

PCX) 

F{z;t)= e- pz c (t + p)dp. (1.25) 
Jo 

It should be noted, however, that formula (1.25) has rather formal meaning. In practice, there 
are situations when direct application of this formula may be problematic, if, e.g. integral in 
rhs (1.25) diverges. As a simple example consider the case when c (t) has the expression 

oo 

c o(t) = ^fcexp(z/ fc t) (1.26) 

fc=— oo 

with some complex constants \ik-,Vk- Then 

c n (t) = d 2^ = £ exp{u k t) 

k=—oo 

and the Stieltjes function F(z; t) has the expression 

F(z;t) = Y.c n {t)z-"- l = £ (L27) 

n=0 k=-oo Z Vk 

i.e. the orthogonality measure in this case is located at the points v k of the complex plane with 
the corresponding concentrated masses M k (t) = fj, k exp(v k t). Formula (1.27) is a special case 
of the well known result for the restricted Toda chain: if dp(x) is an orthogonality measure for 
the orthogonal polynomials P n (z;0) (i.e. for initial value of time t = 0), then for arbitrary t 
the measure will be [4], [18] 

dp(x; t) = const exp(xt) dp{x) (1-28) 

where the constant in rhs of (1.28) is not essential and is needed only to provide the nor- 
malization condition for the measure. Formula (1.27) corresponds to the special case of the 
measure 

oo 

dp{x) = £ p k b[x-v k )dx (1.29) 

k=— oo 
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On the other hand, the generating function t) has the expression 



00 r (t)n n 00 

Hp; t) = E -^f- = E ^ exp(i*(p + *)) (i-so) 



n=0 n - fc=-oo 



If one applies formula (1.23) for the Laplace transform to the function given by (1.30) 
we get 

F(z;t)= E e" fc Vfc / e-"^-"*)dp (1.31) 

Doing "naively" we can put 

/ e- p{z - Vk) dp = (z- v k y x (1.32) 

JO 

in (1.31) and obtain desired formula (1.27) for the Stieltjes function F(z;t). However, formula 
(1.32) is correct only if Re(z) > Re(u k ). For Re(z) < Re(u k ) we may still use the same formula 
(1.31) because it corresponds to the true formal series (1.27). 

Sometimes the following trick will be useful. Consider the modified function co(t) = co(it) 
(i.e. just pass to the "imaginary time"). Construct the Hankel determinants D n (t) form the 
new moments c n (t) = i n c n (it). Clearly, D n (t) = i n ^ n ~^ D n (it). Define corresponding recurrence 
coefficients b n (t),u n (t). We have 

b n (t) = ib n (it), u n (t) = -u n {it) (1.33) 

Obviously they will satisfy the restricted Toda chain equations (1.1). Construct corresponding 
orthogonal polynomials 

P n {z;t)=i n P n (z/i;it) (1.34) 
These "modified" orthogonal polynomials satisfy the recurrence relation 

P n+1 (z; t) + b n (t)P n (z; t) + u n P n ^(z; t) = zP n (z; t) (1.35) 

with the initial conditions 

P (z;t) = l, P 1 (z;t) = z-b (t) 

The Stieltjes function F(z; t) for these polynomials is defined as 

00 00 
F(z;t) = Y,c n (t)z- n - 1 = Y,i n c n (it)z- n - 1 = -iF(z/i,it) 

n=0 n=0 

Assume now that the function c (t) takes real values on the real axis t and is periodic with 
some real period c (t + T) = c (t). Then we can present c (t) (at least formally) in terms of 
the Fourier series 

00 

Co(t)= E fikexp(2mkt/T) (1.36) 

k=— 00 
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In this case the measure for the modified polynomials P n (z; t) is purely discrete and is located 
on the uniform grid 

z k = 2ixik/T, fc = 0,±l,±2,... (1.37) 

on the imaginary axis. Hence, for real and periodic functions co(t) we obtain orthogonal poly- 
nomials P n (z; t) with the measure located on the imaginary axis. 

Returning to the initial orthogonal polynomials P n (z;t) we easily obtain 

oo 

F ^t)= E - _ 2 l k/T eM^kt/T) (1.38) 

fc=— oo / 

We see that in this case the Stieltjes function F(z; t) corresponds to a purely discrete measure 
located on the real axis at points 

z k = 2irk/T, k = 0, ±1, ±2, . . . (1.39) 

with the corresponding discrete masses 

M k (i) = n k exp{2irkt/T) (1.40) 

The measure is well defined if conditions for the moments 

oo oo 

Cj (t)= ]T M k (t)z{= J2 A** exp(2nkt/T)4 < oo (1.41) 

fc=— oo fc=— oo 

hold for all nonnegative integers j (at lest for some values of t). 

The trick with passing to the purely imaginary time t — > it will be especially useful if the 
Fourier series for the modified function co(t) is known. Then we can restore information on the 
initial measure using formula (1.38). 

2. Elliptic functions. 

In this section we recall basic properties of the Weierstrass elliptic functions which will be needed 
in further analysis [2], [26]. The Weierstrass function p(z;g 2 ,g3) depends on the argument z 
and two parameters #2, #3 (the so-called invariants). It satisfies the differential equation 



'2 



(z) = Ap 3 (z) - g 2 p(z) -g 3 = A(p{z) - e 1 )(p(z) - e z ){p{z) - e 3 ) 



where the parameters satisfy the restriction e\ + e 2 + e 3 = 0. The function p(z) is double- 
periodic: 

p(z + 2u) = p(z + 2J) = p(z), 
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where the periods 2u,2u' are assumed to satisfy the condition Im(u'/u) > 0. If #2, #3 are 
known then the periods 2u, 2oj' can be calculated by a standard procedure in terms of the 
elliptic integrals of the first kind [2]. 

It is convenient to introduce the notation [2] 

UJl — UJ, U 3 = u', ijJ 2 = —oj — uj' 

There is a relation between e, and cuf 

p(uj k ) = e k , k = 1,2,3 (2.1) 
The Weierstrass zeta function ((z) is an odd function £(— z) = —((z) defined as 

CO*) = -P(z) 

The function ((z) has simple poles at the points 2mw + 2m'uj', where m,m' are arbitrary 
integers. In contrast to p(z), the zeta function ((z) is quasiperiodic: 

C(z + 2u k ) = C(z) + 2rj k , k = 1,2,3 (2.2) 

where 

Vk = CM 

There are useful relations for r] k : 

Vi + V2 + Vs = (2.3) 

and 

rj 2 uJi - f]\L)2 = in/ 2 

(there are two similar relations which are obtained from the last relation by a cyclic permutation 
of 1,2,3). 

The Weierstrass sigma function a(z) is an odd function cr(— z) = —cr(z) defined as 

It has simple zeroes at the points 2mui + 2m'uj' . The sigma function has quasi-periodic property 

a(z + 2u k ) = - exp(2r] k (z + uj k )) a(z) (2.5) 

There is a simple formula connecting p{z) and a(z) [2]: 

a(u + v)a(u — v) 



a 2 (u)a 2 (v) 



= p(v) - p(u) (2.6) 



Apart from the function a(z) one can define functions <Tk(z), k = 1, 2, 3 by the formulas [2] 

cr a (z) = — + ^ exp(-zr] a ), a = 1,2,3 (2.7) 
a(uj a ) 

The functions a a (z) are convenient when passing from the Weierstrass to the Jacobi elliptic 
functions. Indeed, we have [26] 

sn(u; k) = (ei - e 3 ) 1/2 ^-, cn(u; k) = ^4, dn( M ; fc) = (2-8) 

03(2) 03(2) 03(2) 

where 

u = {e 1 -e 3 )^z, k 2 = 6 -^^ (2.9) 
The parameter fc is called the elliptic modulus. The parameter 



ei - e 2 
ei - e 3 



jf = (1 - A; 2 ) 1 / 2 = 

is called the complementary modulus [26]. The values 

K = y/ei - e 3 ui, K' = iy/e[ - e 3 uj 3 

are complete elliptic integrals of the first kind [26]. 

We need also expressions of the Weierstrass functions for the value z = 0J\/2: 

p(wi/2) = p(3wi/2) = ei + (ei - e 3 )k' , p(uj 3 + wi/2) = ^ ' _^~J 2 , 

p'(wi/2) = -2fc'(l + fc')(ei - e 3 ) 3/2 , p"(cu 1 /2) = 4(ei - e 3 ) (2(d - e 2 ) + 3eifc') (2.10) 

2 p'(a;i/2j 
?7i + Vei - e 3 (fc' + 1), 

((u 3 + wi/2) = r] 3 + rii/2 + v^7^(l - k')/2. (2.11) 

Note that the choice of an appropriate sign in front of square roots in formulas (2.10) and (2.11) 
is not trivial problem and depends on location of the parameters e\, ei, e 3 in the complex domain. 
However in our further analysis we will use the "canonical" choice of these parameters: they 
are real and ordered as e 3 < ei < e±. Then all the square roots are assumed in the arithmetic 
meaning. 

Apart from the Weierstrass zeta function ((z;gi,g 3 ) sometimes the Jacobi Zeta function 
Z(z; k) is more convenient. The relation between these function is [26] 

cw = 3l + vsr^W) + m ("-*>yn , (2.i2) 

io\ { sn{u, k) J 
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where the same relations (2.9) are assumed. 

The Jacobi Zeta function is purely periodic with respect to the period 2K: 

Z(u + 2K) = Z{u) 

and quasi-periodic with respect to the period 2iK': 

Z(u + 2iK') = Z{u) - 

It possesses a remarkable " addition theorem" [26] 

Z(u + v) — Z{u) + Z{v) — k 2 sn(u)sia(v)sn(u + v) (2-13) 

Note the useful relations 

Z(u + K) = Z(u)-k'^P (2.14) 

z(u+K + iK ^ z(u) -^-%. ( , 16) 

The Jacobi Zeta function can be expressed in terms of the incomplete elliptic integral E(u) of 
the second kind [26] 

Z(u) = E(u)-u^ (2.17) 



where 



E(u) = [ U dn 2 (t)dt 
Jo 



and E = E(K). The following formulas allow to express the functions ((z + Ui), i — 1, 2, 3 in 
terms of the Jacobi Zeta function: 



C(z + ui) = rji + — + y/e[ - e 3 < Z(u, k) - 



sn(w, k)da(u, k) 
cn(w, k) 



(2.18) 



C(z + oj 3 ) = 773 + — + VeT r e^Z(u,k), (2.19) 

L0\ 



((z + ^ + lus) = -r/ 2 + — + Vei - e 3 | Z(w, fc) - k 2 

LUl { 



sn(u, k)cn(u, k) 
dn(w, k) 



(2.20) 
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3. Toda chain solution and the corresponding orthogonal 
polynomials 

We present here an explicit solution of the restricted Toda chain 
Lemma 1 Put 

u n (t) = w 2 n 2 (p(w(t + (3)) - p(nw(t + (3) + q)) (3.1) 

and 

b n (t) =m + w(n+ l)C(w(n + l)(t + /3) + q)- wn((wn(t + /3)+q)- (2n + l)w£(w(t + /?)), (3.2) 

where w, (3, q, /ii are arbitrary complex parameters. We will also assume that cui, u 3 are arbitrary 
independent periods corresponding to the arbitrary parameters e±, e 2 , e 3 with the only condition 
ei + e 2 + e 3 = 0. 

Then u n (t),b n (t) satisfy the restricted Toda chain equations (1.1) 

Proof. In order to verify the first equation in (1.1) we present u n (t) in an equivalent form 

22 o-((n + l)w(t + (3) + q)a((n - l)w(t + g) + q) 
Un{t) ~ WH o*(nw(t + /3) + q)a\w{t + /?)) ^ 

using formula (2.6). Then by (2.4) the expression u n /u n can be presented as a sum of the 
Weierstrass zeta functions and we arrive at the first equation of (1.1). The second equation 

(1.1) is satisfied by the formula p(z) = —('(z). 

It is directly verified that the recurrence coefficients u n (t),b n (t) are double-periodic with 
the periods 2u/w,2u' /w: 

u n (t + 2u/w) = u n {t + 2u'/w) = u n (t), b n (t + 2u/w) = b n (t + 2J jw) = b n {t) 

(periodicity property for u n (t) is obvious and periodicity for the coefficients b n (t) follows from 

(2.2) and (2.3)). Thus both u n (t) and b n (t) are elliptic functions in the argument t. 

Using formulas (1.19) and (2.4) we can restore the function c (t). It is easy to verify that 

. a(w(t + P) + q) , . , „. 

Co{t) = a(q)a(w(t + (3)) eXP ^ t + ^ (3 ' 4) 

where /io is an arbitrary constant. 

We thus obtained some new family of orthogonal polynomials P n (z; t) which can be defined 
through given recurrence coefficients u n (t),b n (t). The Stieltjes function F(z,t) (and hence, in 
principle) the orthogonality measure for these polynomials can also be found explicitly from 
formula (1.25) because the function c (t) is given explicitly by (3.4). 
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The obtained orthogonal polynomials P n (z; t) contain several free parameters (say w, q, (3, fii, t 
and elliptic parameters g2,9z)- We would like to investigate some simple special choice of these 
parameters when our polynomials are a generalization of already known families. Note that 
the parameter (3 is inessential: it describes a shift of the argument t — > t + (3. Nevertheless, 
we will keep this parameter for convenience, assuming that the argument t takes real or pure 
imaginary values. We will assume also that q ^ 0- Indeed, the case q = corresponds to 
some degeneration: b (t) = const, Ui(t) = and c (t) becomes a pure exponential function 
c = exp(/xt) in this limit. 

In what follows we put q = Uj, (3 = cu k /w, where k,j are arbitrary noncoinciding integers 
from the set 1, 2, 3. We denote also uj\ = —ujj — uj k . Then using (quasi)periodicity properties of 
the Weierstrass functions p(z)X( z ) we find the expression for u n {t): 

u 2n(t) = 4:W 2 n 2 (p(wt + Uk) — p(2wnt + ujj)) 

u 2n +i(t) = w 2 (2n + l) 2 (p(wt + oo k ) - p(w(2n + l)t + w,)) (3-5) 

and for the coefficients b n (t): 

b 2n (t) = /xi + w {(2n + l)C((2n + l)wt - uj t ) - 2n((2nwt + Uj) - (4n + 1)CM + cu k ) + 2nr] k } , 

= A*i + w {{2{n + l)C((2(n + l)wt + Uj) - (2n + l)C((2n + l)wt - u{)- 
(4n + 3)((wt + tu k ) + (6n + 4)^} (3.6) 

We will assume also that fi — and 

/ii = -wrjj, (3.7) 

Indeed, the parameter /i is inessential and can be chosen arbitrary whereas the parameter /ii 
leads only to a trivial shift of the recurrence coefficient b n , hence we can put /xi to a prescribed 
value without loss of generality. 

In order to find the orthogonality measure for the obtained polynomials P n (z;t) we will 
assume that the parameter w is real. Among all 6 possible choices q = Uj, (3 = uj k /w of the 
parameters q, (3 we consider only the two cases: 

(i) if (3 = cui/w, q = cu 2 then form (2.7), (3.7) and (2.8) we find that 

c (t) = Ci/cn(w v /e i — e 3^ k) (3.8) 

where 
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(ii) if (5 — cui/w, q = ojz then quite analogously we find 

c (t) = C 2 dn , (w ^ 61 ~ 6 f § = C 2 dc( Wv ^f; fe) (3.9) 
cnfw^/ei — e 3 t; fc) 

where 

Note that the constant factors Ci,C2 are in fact inessential (the orthogonal polynomials 
P n (z;t) as well as recurrence coefficients b n (t),u n (t) do not depend on these constants) and we 
can put C2 = C\ = 1 . We thus have that Co(i) = l/cn(wty/ei — e 3 ; k) for the case (i) and 
c (t) = dc(w v / ei — e 3 t; fc)) for the case (ii). 

We will restrict ourselves with the case when all are real and distinct, say e\ > e 2 > e 3 
(this is the usual convention [2]). Then it is well known [2] that for the real values of z the 
functions p(z),£(z),cr(z) take real values. For purely imaginary values of z the function p(z) 
takes real values whereas functions ((z),a(z) take purely imaginary values. The period 2a>i is 
real and the period 2a> 3 is purely imaginary. This means that the fundamental parallelogram 
of the elliptic functions in this case is a rectangle [2]. We have also that both k and k' are real 
parameters taking values in the "canonical" interval < k, k' < 1. Hence all values of Jacobi 
elliptic functions sn(x), ca(x), da(x) are real for real x. The functions cn(x),da(x) take also 
real values for purely imaginary values of x. 

Now we are ready to calculate the orthogonality measure for the cases (i) and (ii). 



4. The orthogonality measure and recurrence coefficients 
for the case (i) 

According to considerations of the first section, introduce the new function co(t) = c (it). 
For the first case (i) we have 

c (t) = l/cn(iwy/ei — e 3 i; k) = cn(wy/ei — e 3 t] k') (4.1) 

The function cn(w v / ei — e 3 t; k') is real on the real axis and periodic with the period T = 
— — . The Fourier series for this function is well known [261 

TT °° 1 

cn(«Vei - e 3 t; k') = — Jl ^i/ 2 + v i/2-n e M™(™ ~ lfflwy/el - e 3 t/K'), (4.2) 



where 

v = exp(—7iK/K') 
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Hence the polynomials P n (z;t) have a purely discrete orthogonality measure located at the 
points 

x n = y (2n - 1) = 63 (2n - 1), n = 0, ±1, ±2, . . . (4.3) 

with the corresponding concentrated masses 

n[) ~k'K' v n-l/2 + v l/2-n ^> 

Thus orthogonality relation for the polynomials P n (z; t) looks as follows 

oo 

M s {t)P n {x s -t)P m {x s -t) = h n {t) 5 nm (4.5) 

s=— oo 

It is interesting to determine conditions under which the measure is well defined, i.e. that all 
the moments are finite 

oo 

Cj(t) = J2 M s(tH < oo, j = 0,1,2,... (4.6) 

s=— oo 

From the explicit expression (4.4) it is easily seen that condition (4.6) will hold provided 

K <t< K (4.7) 



uV e i - e 3 w^e x - e 3 

When the parameter t belongs to this interval all the moments are well defined. If t — > ± — — 
then co(t) — > oo as is easily seen from (4.1). Hence when t approaches the endpoints of the 
interval (4.7), the moments c n (t) tend to infinity and the measure becomes not well defined. 
From (4.4) it is clear that for all values of t from the admissible interval (4.7) the concentrated 

masses are positive M s (t) > 0, s = 0, ±1, ±2, This means that we indeed deal with a 

positively defined purely discrete measure on the whole real axis. 
Moreover it is easy verified that 

£ M s (t)xl = Cj (t) = — c (t) (4.8) 

s=-oo at 

where Co(t) = l/cn(wty/ei — e?,\ k). Indeed, formula (4.8) follows directly from the Fourier 
series (4.2) by j-fold differentiation with respect to t. Formula (4.8) shows that the obtained 
measure is "true", i.e. it gives the prescribed moments Cj(t) for all j — 0, 1, . . .. 

Consider the recurrence coefficients b n (t),u n (t) for the orthogonal polynomials P n (z;t) cor- 
responding to the function c (t) defined by (4.1). 

It is clear that both b n (t) and u n (t) are real for all t from the admissible interval (4.7). 
Indeed, the Hankel determinants D n (t) are real because these are constructed from the matrix 
with real entries 

a-ii = dt i +j » = 0,1,..., n- 1 
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Hence the normalization coefficients h n (t) are real as well. The same is true for coefficients 
b n (t),u n (t) obtained from h n (t) by formulas 

b n (t) = j i , u n = h n /h n _ 1 
Explicitly we have for the recurrence coefficients b n (t) (see (3.6)) 

M*) =w{(2n + l)C((2n + l)wt - uj 3 ) - 2n((2nwt + u 2 ) - {An + l)({wt + + 2nn x - r? 2 } , 

= ™ {(2(n + l)C((2(n + l)wf + u 2 ) - (2n + l)C((2n + l)wt - u 3 )- 
(4n + 3)({wt + wi) + (6ra + 4)771 - r? 2 } (4.9) 

or, in terms of the elliptic Jacobi functions 



b 2n {t) = Wei - e 3 {(2n + l)Z((2n + l)u) - 2nZ(2nu) - (4n + l)Z(u)- 

sn(-u)dn(w) 1 



2n£; 



2 cn(2rax)sn(2mi) 
dn(2rax) 



+ (4n + 1) 



emu) 



(4.10) 



uVd - e 3 I (2n + l)Z((2n + l)u) - 2nZ(2nw + X) - (4n + 1) + AT + iAT') + J 



= ^Vei - e 3 {2(n + l)Z(2(n + l)u) - (2n + l)Z((2ra + l)u) - (4n + 3)Z(u)- 



2 cn(2(n + l)w)sn(2(n + i i ^ sn(-u)dn(w) 



2(n + l)£: 2 vv ' // \ v - ' ' + (An + 3) > = (4.11) 

v ; dn(2(n + l)u) V ; cn(u) J V ; 

Wei - e 3 |2(n + l)Z(2(n + l)u + AT) - (2n + l)Z((2n + l)u) - {An + 3) (z(u) + } 
where 

u = wt y/e[ — e 3 

For the recurrence coefficients w n (t) we have expressions in terms of the elliptic Jacobi 
functions 



, ., , , 1 | fc2 sn 2 (2rm) \ 

cn 2 (w) dn 2 (2n«) / 



«2n(t) = An 2 w 2 (e 1 - e 2 ) — + k 



u 2n+l {t) = (2n + 1) V(ei - es ) (V + dn 2 ((2n + 1)«)^ (4.12) 

It is seen from (4.12) that for all t from the admissible interval the recurrence coefficients 
u n [t) are bounded and strictly positive u n (t) > 0. As is well known from general theory 
of orthogonal polynomials [8] the property u n > for all n > (together with reality of 
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the coefficients b n ) guarantees existence of a positive measure on the real axis. We already 
constructed this measure explicitly (4.5). 

The remaining question is about uniqueness of the moment problem for the case (i). Indeed, 
we have constructed explicitly the orthogonality measure (4.5) corresponding to the moments 



cn(t) = — — , n = 0,l,2,..., (4.13) 

dt n [ cn^wtvei — e 3 ; k) J 

where t is assumed to belong to the admissible interval (4.7). 

But in principle, it is possible that this measure is not unique. Such situation is known 
as indeterminate moment problem [1], [20]. In more details this means the following. Assume 
that real moments c n are given and all corresponding the Hankel determinants D n > are 
positive. This condition is equivalent to positivity of the recurrence coefficients u n > for 
n — 1,2, .. . and in turn, it guarantees existence of a positive orthogonality measure on the 
real line — oo < x < oo (so-called the Hamburger moment problem [20], [9]). If this measure 
is unique (up to a normalization condition) then the Hamburger moment problem is called the 
determinate. If there exist at least two different orthogonality measures then the Hamburger 
moment problem is called indeterminate. In case of the indeterminate Hamburger problem 
there exists infinitely many different measures (see [1], [20] for details). 

Finding criteria for determinacy of the Hamburger moment problem is a nontrivial problem 
[9]. However, there is a simple sufficient condition proposed by Carleman [20]: if 

oo 

£< 1/2 = oc (4.14) 

n=l 

then the Hamburger problem is determinate. Of course, in (4.14) the arithmetic value of the 
square root ^fu^ is assumed. 

We now show that the Hamburger problem for the moment problem (4.13) is determinate. 

From (4.12) it follows that for all n and for fixed t from the admissible interval we have the 
inequalities 

< U 2n (t) < A(t) 

where 

1 k 2 dn 2 (w) 



cn 2 (w) k' 2 k' 2 cn 2 (u) 
is a fixed positive parameter (depending on t but not on n) . Hence we have 



E 

n=l 



-1/2 

u 2n = oo 
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diverges. From the similar considerations it follows that 

oo 

E-l/2 
U 2n +1 = °o 

n=l 

Hence the Carleman condition (4.14) holds and we indeed have the determinate moment prob- 
lem. This means that the discrete measure (4.5) is the only providing orthogonality of the 
polynomials P n (z; t) on the real axis. 

When t approaches the endpoints of the admissible interval (4.7) the recurrence coefficients 
b n (t),u n (t) tend to infinity. This "explosion" of the recurrence coefficients explains correspond- 
ing "explosion" of the orthogonality measure when t tends to he endpoints of the admissible 
interval. 

Put now t — (the middlepoint of the admissible interval) . Then it is seen (due to property 
f]\ + f}i + f]z = 0) that the recurrence coefficient b n (0) vanishes 

6 n (0) = 0, n = 0,l,2,... (4.15) 

and we arrive at a class of so-called symmetric orthogonal polynomials with the recurrence 
relation [8] 

P n +i(z) + u n (0)P n ^(z) = zP n {z) (4.16) 
For the recurrence coefficients u n (0) we have 

u 2n (0) = 4w 2 n 2 (e 1 -e 2 ), u 2n+1 (0) = W 2 (2n+l) 2 ( ei -e 3 ). (4.17) 

These recurrence coefficients correspond to the orthogonal polynomials introduced by Carlitz 
[7]. In turn, orthogonality relation for these OP follows from the remarkable result by Stieltjes 
on presenting of the Laplace transform of the Jacobi elliptic functions in terms of continued 
fraction (for modern treatment of this and related examples see e.g. [16]. Today the orthogonal 
polynomials introduced by Carlitz are called the Stieltjes-Carlitz orthogonal polynomials related 
with elliptic functions [8]. For further details concerning these polynomials see [15]. 

The orthogonality measure for polynomials corresponding to the Stieltjes-Carlitz case (4.17) 
is obtained from our measure by putting t — 0. We thus have the orthogonality relation 

oo 

J2 M s (0)P n (x s ; 0)P rn (x s ; 0) = ^(0) t ( 4 - 18 ) 

s=— oo 

where the support of the measure is the same, i.e. the points x s have the same expression (4.3) 
and the concentrated masses are 
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This measure was discovered by Stieltjes (Carlitz showed that this measure provides orthog- 
onality of the corresponding polynomials P n (z)). For the Stieltjes- Carlitz polynomials the 
Hamburger moment problem is obviously determined because t — belongs to the admissible 
interval. 

We see that the Stieltjes-Carlitz polynomials appear naturally as a very special case of the 
"elliptic Toda polynomials" corresponding to "zero time" condition t — 0. 
We can finally summarize all these results as the 

Theorem 1 Assume that e 3 < e 2 < t\ are arbitrary real parameters with the condition e x + e 2 + 
e 3 = 0. Assume that the recurrence coefficients are given by formulas (4-9), (4-12) with arbitrary 
positive parameter w. Assume also that the parameter t belongs to the admissible interval (4-7). 
Then the corresponding orthogonal polynomials P n (x; t) are orthogonal on the uniform grid (4-3) 
on the real axis with the concentrated masses given by (4-4)- The corresponding moments c n (t) 
are given by (4-13). The moment problem is determinate. 

5. The orthogonality measure and recurrence coefficients 
for the case (ii) 

Consider now the case (ii). We have analogously 

Co(t) = dc(wu-\/ e i — e 3 t; k) = dn(w v /e i — e 3^; k') (5.1) 
The Fourier series is well known [26] 

71 00 1 

dn(w y /e 1 - e 3 t; k') = — ^ - n exp(7rinwt y /e 1 - e 3 / K') (5.2) 

1l = — oo ^ I V 



with the same expression for h as for the case (i). From considerations of the first section we 
see that corresponding orthogonal polynomials P n [z\ t) have purely discrete measure located at 
the points 



K< 

Corresponding concentrated masses are 



-, n = 0,±l,±2,... (5.3) 



2tt 

M n(t) = 777? n T - n \ ^M^n^^t/K') (5.4) 
K [y + v ) 



The admissible interval for t is the same as for the case (i): 

K <t< K 
wyfei - e 3 w^ei - e 3 
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As for the case (i) it is easily verified that inside the admissible interval the recurrence coeffi- 
cients b n (t),u n (t) are real and u n > which guarantees positivity of the measure: M n (t) > 
for all n — 0, ±1, ±2 and for all t from the admissible interval. 
The coefficients u n (t) has the expression 

u 2n (t) = Wn\e x - e 3 ) (k' 2 + dn 2 (2nu)) 

u 2n+l (t) = w 2 (2n + l) 2 ( ei - e 2 ) [-^ + k 2 E^LtM] (5.5) 

\cn 2 (u) dn ((2n + l)u) ) 

The recurrence coefficients b n (t) are expressed as 

hn(t) = wVei - e 3 {(2n + l)Z((2ra + l)u + 
2nZ(2nw) - (4n + l)(Z(u + AT + iK') - m/{2K))} 

(5.6) 

Wi(*) = ^v^i^i {2(n + l)Z(2(n + l)u)- 

(2n + l)Z((2n + 1)« + AT) - (4n + 3)(Z(u + K + - m/(2K))} 

From (5.5) it is clear that the coefficients u n (t) are strictly positive u n (t) > for any fixed 
value of the parameter t from the admissible interval. 

From the same considerations it follows that the Hamburger moment problem for the mo- 
ments 

r . m _ dn j dn(w yei - e 3 t; fc) 1 _ 
nl J ~ dt" \ cniwty^^t; k) } ' n " U ' X ' ^' • • • 

is determinate for any value of the parameter t from the admissible interval (4.7). 

When t — then again, as in the case (i) the diagonal recurrence coefficients are zero 
6 n (0) = and 

u 2n = 4u?V(ei - e 3 ), w 2 n+i = w 2 (2n + l) 2 (ei - e 2 ) (5.7) 

Note that the recurrence coefficients (5.7) are obtained from the corresponding coefficients 
(4.17) of the case (i) by a simple transposition e 2 <-> e 3 . These recurrence coefficients (5.7) 
correspond to the second class of the Stieltjes-Carlitz orthogonal polynomials [7], [15] arising 
from the Laplace transformation of the elliptic function dn(i; k'). 

We see, that just as in the case (i), the special choice of the parameter t — leads to already 
known Stieltjes-Carlitz orthogonal polynomials. For arbitrary values of the parameter t from 
the admissible interval we obtain new orthogonal polynomials with explicitly known recurrence 
coefficients and positive discrete measure on a whole real axis. 

We considered only two possible choices of the parameters (3, q (f3 — Wi/w, q = oo 2 and 
j3 = wi/w, q = u>s) because for these two choices we get polynomials P n (z;t) having positive 
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orthogonality measure on the real axis. It seems that all other values of the parameters (3, q 
(for real values w) do not lead to polynomials with a positive measure. Nevertheless, as we will 
see in the next sections, in the degenerate cases of elliptic functions there are more possibilities 
for these parameters when the measure appears to be a positive. 



6. Special choices of the parameter t 

The expressions (3.5) and (3.6) show that for generic value of t the recurrence coefficients 
b n (t),u n (t) are transcendent functions in n. We already saw that for special choice t = we 
obtain b n (0) = and u 2n (0) as well u 2n+ i(0) are quadratic polynomials in n. 
Put now 

where M < N are mutually prime positive integers. Present the number n in the form 

n = Nr + s, r = 0, 1, . . . , s = 0, 1, . . . , N - 1 

Then for fixed s = 0,1,..., N — 1 the recurrence coefficients b 2n (t), b 2n+ i(t) will be linear 
function in r and the coefficients u 2n (t), u 2n+ i(t) will be quadratic polynomials in r . 

Indeed, from formulas (3.5), (3.6), using periodicity of the Weierstrass functions, we get 

u 2Nj+ 2s = 4w\Nr + s) 2 (e - e 1 {s)) , (6.2) 
u 2N r+2s+i = w 2 (2Nr + 2s+ l) 2 (e - e 2 (s)) , 

where 

e = p((M + N)ujN), ei(s) = p(uj + 25^^^!/^), e 2 (s) = p(u t + (2s + tfMujN) 

In these formulas j — 2, 1 — 3 for the case (i) and j = 3, / = 2 for the case (ii). 
Analogously, for the coefficients b 2n (t) and b 2n +\{t) we obtain 

b 2N r+2s = w{-r}j + 2i ll {Mr + Nr + s)) + w(K 1 (s)(2Nr + 2s + 1)) - 
2K 2 (s)(Nr + s) - K (A(Nr + s) + 1)) (6.3) 

b 2N r+2s+i = w(-r}j + 2r} x {Mr + 3Nr + 3s + 2)) + w(K 3 (s)(2Nr + 2s + 2)) - 
K 4 (s)(2Nr + 2s + 1) - K (A(Nr + s) + 3)) (6.4) 

where 

k = C((M + N)ui/N), «i(s) = C(-wj + (2s + ^MujN), k 2 (s) = ((ujj + 2sMcj 1 /N) 
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k 3 (s) = ((uj + 2sMu)i/N), Ki{s) = ((-cui + sMui/N) 

We see that indeed u 2 n t+2s , u 2 N r+2s+ i are quadratic in r and b 2 ^ r+2si b 2 Nr+2s+i are linear in n. 

The constants 6j(s), can be found in a less or more "explicit" form only for several values 
of N, M. We already considered the case t — which corresponds to M = 0, N — 1, s — 0. 
Another simplest case corresponds to the choice 

M = 1,N = 2 

i.e. we choose t — In this case the parameter s can take only 2 values s = 0, 1. 
The zero moment becomes now 



co(t) = l/cn(X/2; fc) = 
Using relations and (2.11) we can calculate the recurrence coefficients b n (ui/(2w)) 



b 4n (uj 1 /(2w)) = w Vei - e 3 (2n(k' + 3) + 1) (6.5) 

b 4n+1 (uj 1 /(2w)) = w Ve~i -e 3 (2n(3k' + 1) + 2k' + 1) (6.6) 

b An+2 (oj 1 /(2w)) = w Vei -e 3 (2n(3fc' + 1) + 4A;' + 1) (6.7) 

&4n+3(^i/(2u0) = w Vei - e 3 (2n(k' + 3) + 2A/ + 5) (6.8) 
Similarly using relations (2.10) we obtain expressions for the coefficients u n (ui/ (2w)); 

u 2n+1 {u 1 /{2w)) = w 2 (ei - e 3 ){2n + l) 2 2k' (6.9) 

u 4n (wi/(2iu)) = 16w 2 ( ei - e 3 )n 2 k'(l + k') (6.10) 

u An+2 {uj 1 /{2w)) = Aw 2 {e 1 - e 3 )(2n + 1) 2 (1 + k') (6.11) 



Note that the combination W\/e\ — e 3 plays the role of a scaling parameter, hence we can 
put w — 1, ei — e 3 = 1 without loss of generality. 
We then have 



ll + k' 



21 



Theorem 2 Assume that < k' < 1 is an arbitrary parameter. Let the recurrence coefficients 
b n be defined as 

(k! + 3)n/2 + 1, i/ n = (mod 4); 

(3&/ + l)n/2 + (A;' + l)/2, i/ n = 1 (mod 4); 

(3A/ + l)n/2 + fc', if n = 2 (mod 4); 

(Jfe' + 3)n/2 + (A;' + l)/2 if n = 3 (mod 4) (6.12) 

and £/ie recurrence coefficients u n be defined as 

k'(k' + l)n 2 , i/ n = (mod 4); 
2A;'n 2 , i/ n = l,3 (mod 4); 

(1 + A;> 2 , if n = 2 (mod 4) (6.13) 

T/ien corresponding monic orthogonal polynomials P n (x) are orthogonal with purely discrete 
measure on the real line 

oo 

^ ] M s P n (x s ) P m (x s) h n S nm , (6.14) 

s=— 00 

where the grid of orthogonality is 

n(s - 1/2) 
Xs ~ K> 
and corresponding concentrated masses are 

n v (s-l/2)/2 

M ' = UK' *-W + v = zM-*K/K>) 

The normalization constants h n are 

h n = c Mi-u 2 •■•«« 

where 

c = 



'1 + Jfe 



k' 



7. Degenerated cases 

Consider degenerated cases of obtained orthogonal polynomials. These degenerated cases arise 
when two or three of the parameters t{ coincide. Using our choice e± > e<i > e 3 we see that 
there are two possibilities when two of the parameters coincide: 
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(i) e\ = C2 = a, = —2a, where a is a positive parameter. In this case the real period 2uj\ 
tends to infinity, whereas the imaginary period remains finite 



Til 

Without loss of generality we can take a = 1/3 (changing of a leads only to scaling of the argu- 
ment z of corresponding functions). The Weierstrass functions are then reduced to hyperbolic 
ones, e.g. 

p(z) ^1/3+—^— 
sinn (z) 

The modular parameter becomes k = 1 and the Jacobi elliptic functions become hyperbolic as 
well: sn(z; 1) = tanh(z), cn(z; 1) = da(z; 1) = 1/ cosh(z). The imaginary period in this case is 
2a> 3 = in. 

(ii) C2 = e3 = —a, e\ = 2a with some positive parameter a. Then the imaginary period 2u!s 
becomes infinity whereas the real period is finite 
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Again we can put a = 1/3, then the Weierstrass function p{z) is degenerated to trigonometric 
form: 

1 

p(z) -> -1/3+ , 2 

sin (z) 

The modular parameter k = in this limit and we have sn(^; 0) = sm(z), cn(z; 0) = cos(z), dn(z; 
1. 

Consider first the degenerated cases (i) and (ii) of the elliptic polynomials obtained in 
the previous section. In the "hyperbolic" limit k — 1 the functions c (t) becomes the same 
c (t) = cosh(t). This case is non-interesting because it corresponds to degeneration of the 
orthogonal polynomials: the Hankel determinants become zero D n (t) for infinitely many n. In 
turn, this corresponds to zero recurrence coefficients for even 2n: U2 n (t) = as can be seen 
from explicit formulas for u n (t) in the hyperbolic limit. 

In the trigonometric limit k — again both functions c (t) coincide: c (t) = 1/ cos(t). 
This function c (t) corresponds to some elementary solutions of the restricted Toda chain; 
corresponding orthogonal polynomials P n (z; t) coincide with a special case of the Meixner- 
Pollaczeck polynomials. 

Indeed, for this case it is elementary verified that the recurrence coefficients have the ex- 
pression: 

u n(t) = — 27, Kit) = (2n + 1) tan* (7.1) 
cos z t 
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On the other hand the monic Meixner-Pollaczeck polynomials [13] 

(here (a) n = a(a + 1) ... (a + n— 1) is the standard shifted factorial) depend on two parameters 
A, and have the recurrence coefficients 

njn + 2A-1) n + \ 

U n = . ■ 2 i > & ™ = 7 T (7-3) 

4 sin tan0 

Comparing the recurrence coefficients we see that our polynomials coincide with (rescaled) 
Meixner-Pollaczeck polynomials with the parameters A = 1/2, = t + 7r/2. 
The explicit expression for our polynomials P n (x;t) is 



P„ (l;() .!^! 2Fl f-«.V2 + ix/2 _ 2 „ 

COS t 



+ e- 2it (7.4) 



From the standard formulas for the Meixner-Pollaczeck polynomials [13] we find that that our 
polynomials P n (x;t) are orthogonal on whole real axis: 

roo 

/ P n {x; t)P m {x; t)W(x; t)dx = h n (t) S nm , (7.5) 

J x=—oo 

where the weight function W(x,t) is 

tx 

W(x,t) = - — ^ — — (7.6) 
v ' ; 2cosh(7rx/2) v ; 

The weight is well defined provided that t belongs to the admissible interval — 7r/2 < t < tt/2. 
When t is inside this interval all moments c n (t) exist and we have the normalization condition 

tx 



W(x; t)dx = I — -rrzdx = c (t) 

-oo y ' ; 7-00 2 cosh vrx/2 v ; 



oo 2 cosh(7ra;/2) ' ' ' cost 

It is instructive to see how the continuous orthogonality relation (7.5) arises in the limiting case 
of the discrete-type orthogonality relation (4.5). Indeed, in the trigonometric limit k — > we 
have K(k) — > ir/2 and K'(k) — > oo. So (recall that we assume e\ = 2/3, e 2 = e 3 = —1/3) from 
(4.3) we see that the grid step Ax(s) = x s+ \ — x s becomes infinitely small and the sum in rhs 
of (4.5) becomes an integral (7.5) after appropriate definition of the continuous variable x. 

Note also that the special case considered in the previous section (recurrence coefficients 
given by (6.12), (6.13)) in the limit k — corresponds to the formulas (7.1) for t = n/A: 

u n = 2n 2 , b n = 2n + 1. 
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Consider now more general class of degenerated solutions corresponding to the case when 
(3 = 0, Hi = wq/3, (J>o — and q is an arbitrary real parameter. We then have in the hyperbolic 
limit (k = 1) 

smh(wt + q) 1 2 

0{ 1 sinh(g) sinh(wt) 2(e 2 " - 1) ^ 1 - e~ 2wt K ' 



Introduce the function 



wt 



C ° 0)(t) "l-e-^-sinh( W t)- (7 ' 8) 

The function (7.7) differs from c °\t) only by adding of a term e~ q / sinh(g) not depending on 
t. We have 

Cn(t) = <$Ht) = *^, n = l,2,... 

Thus all moments corresponding to functions c °\t) and c (t) coincide apart from zero moments. 
Introduce linear functional a it) and a^\t) by their moments 

(a(t),x n ) = c n (t), (^°\t),x n ) = c^(t), n = 0,1,2,... 

We see that the functionals a(t) and a^(t) are related as 

<t(°) (t) = (7(0 + e"V sinh(g)<5 , (7.9) 

where <5 is the Dirac delta-functional, corresponding to inserting a unit concentrated mass to 
the point x = 0: 

Assume that p^(x) is the orthogonality measure for the polynomials P^(z;t) corresponding 
to the functional a^(t): 

/CO 
P^(x;t)P^(x;W°\x;t) = 5 nm 
-oo 

Then the orthogonality measure corresponding to he polynomials P n (x; t) is 

p(x; t) = p (0) (x; t) + e~ q / sinh(g) 5{x) 

Thus indeed the weight of orthogonality for the polynomials P n (x;t) is obtained from the 
corresponding orthogonality weight for the polynomials P^ (x; t) by adding of a concentrated 
mass coth(g) at the point x — 0. 

Now we show that the function c °\t) given by (7.8) generates a special class of the Meixner 
polynomials. Indeed, it is elementary verified that two sequences 

, /N 2w(n + 1 + ne 2wt ) Aw 2 n 2 e 2wt n ^ n 

b n (t) = l , u n = {e2wt _ 1)2 , n = 0,l,2,... (7.10) 
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satisfy the restricted Toda chain equations (1.1) together with the condition bo(t) = (t) / (^\t) . 
Thus the recurrence coefficients (7.10) correspond to orthogonal polynomials P^(z;t) having 
the moments cffl(t) = d n c^\t)/dt n . On the other hand, we can easily identify the recurrence 
coefficients (7.10) with the those for the special class of the Meixner polynomials. 

Indeed, the Meixner polynomials P n (x;j3;c) have two real parameters (3, c and have the 
recurrence coefficients [13] 

n+(n + (3)c c . n iN 

b n = 7 — , u n = - - n (n + P-l) 7.11 

1 — c (1 — c) z 

Explicitly the Meixner polynomials are expressed in terms of the Gauss hypergeometric function 
[13] 

p ; 1 - 1/cJ (7.12) 

where appropriate normalization factor to provide monicity of the polynomials P n {x; (3; c). 

The Meixner polynomials are orthogonal on the uniform semi-infinite grid [13]: 

£ C ^TT P ^ & c ) p ^ & c ) = K 5 ™ ( 7 - 13 ) 

Obviously we should have < c < 1 to provide positivity property of the measure. 

Comparing the recurrence coefficients (7.10) with (7.11) we see that (3 = 1, c— e~ 2wt and 
polynomials P^ \x; t) coincide with the corresponding rescaled Meixner polynomials: 

Pi°\x;t) = (-2w) n P n (-x/(2w);l;e- 2wt ) (7.14) 

Orthogonality relation for the polynomials P^ (x; t) looks as 

oo 

Y,e- 2swt P£\-2ws-t)P£\-2ws-t) = hSHt) 5 nm (7.15) 

s=0 

To provide positivity of the measure for t > we should have w > 0. Thus polynomials 
P^(x; t) are orthogonal on the uniform grid of the negative real axis. 

Return to the polynomials P n (z;t) corresponding to the function (7.8). The recurrence 
coefficients for the polynomials P n (z;t) are obtained from the recurrence coefficients (3.2), 
(3.1) by the limiting procedure e 2 — > e\. 

b n (t) = w(n + 1) coth(w(n + l)t + q) — wn coth(wnt + q) — w{2n + 1) coth(wt) (7.16) 

and 

/ . \ 2 2 sinh(w(n + l)t + q)wah.[w[n- l)t + q) 

u n \t) = w n -07 — r (7.17) 

sinh (wnt + q) sinh (wt) 
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As we already showed, the orthogonality relation for the polynomials P n (x; t) corresponding 
to the function (7.8) is obtained from (7.15) by adding of a concentrated mass at the point 
x = 0. Explicitly we have 

oo 

£ 2e' 2wst P n (-2ws, t)P m (-2ws, t) + MP„(0; t)P m (0; t) = h n 5 nm (7.18) 

where the value of the mass inserted at x = is 

M 



sinh(g) 

It is assumed that w,t > in order to provide convergence of series in lhs of (7.18). 

The normalization constants h n (t) are expressed through the recurrence coefficients (7.17) 

as 

h n (t) = c (t) Ul (t)u 2 (t) . . . u n (t) (7.19) 

Note that for n = m = formula (7.18) gives an obvious identity 

y e -2 W ts + e~ 9 = sinh(wt + q) = 
~^ sinh(g) sinh(g) sinh(wt) 

When q — > oo we see that M — > and c (t) — > si ^( w t) = c o°' ) (^) hence in this limit the 
polynomials P n (z;t) become the ordinary Meixner polynomials (7.14). 

The polynomials obtained by an adding of a concentrated mass at the point x = of the 
orthogonality measure for the Meixner polynomials are called the modified Meixner polynomials 
and were proposed by R.Askey as an interesting object for further investigations [5]. Properties 
of these polynomials were intensively studied in [3] and [6]. In particular it was shown that 
these polynomials satisfy difference equations of finite and infinite order. 



8. Trigonometric limit 

Consider the trigonometric limit when e\ = 2/3, e 2 = = —1/3. Put /3 — 0, //o = 0, \i\ — 
—wq/3. The recurrence coefficients take the form 

b n (t) — w(n + 1) cot(w(n + l)t + q) — wn cot(wnt + q) — (2n + l)w cot(wt) (8.1) 

and 

U n {t) 

and the function c (t) is 



n o sin((n + l)wt + q) sin((n — l)wt + q) 
w n 



Co(t) 



sin (nwt + q) sin (wt) 
sin(wt + q) 



sin(g) sin(wt) 



(8.2) 
(8.3) 
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From the expression (8.2) we see that for real value t,w the coefficient u n (t) cannot have the 
same sign for all n. Hence, corresponding measure for orthogonal polynomials is not positive 
definite. 

Nevertheless, there is one interesting special case leading to a positive definite measure on 
a finite set of points on the real axis. Indeed, put w — 1, q — 7r/2, then 

Co(t) = cot(f) (8.4) 

and 



b n (t) = -(n+ l)tan((n+ l)t) +ntan(nt) - (2n + 1) cot(f), 
u n {t) = n 2 



2 cos((n + l)t) cos((n — l)t) 



cos 2 (nt) sin (t) 
Assume that 

7T 

t = T 



2(N + 2) 

for some positive integer N = 2, 3, Then it is seen that w fc > for 1 < k < N and un+i = 0. 

This condition guarantees that the finite set of polynomials Po(x; r), P\{x\ r), . . . , P/v(x; r) will 
be orthogonal on the set of (real) zeros x s of the polynomial P/v +1 (x; t): 



N 



J2 PsP n (x s ; T)P m (x s ; r) = h n 5 nm (8.5) 



s=0 



with postive weights p s . 

However explicit expressions for zeros x s and the weights p s in this case are still unknown. 



9. Completely degenerated case. The Krall-Laguerre 
polynomials 

Finally, consider the case when all roots coincide e\ = e 2 = e 3 = 0. Then the elliptic functions 
are degenerated to simple rational ones: p{z) = 1/z 2 , ((z) — 1/z, a(z) = z. 
Assuming again w — 1, (5 — we have 

c (t) = ^ = l/t + l/g (9.1) 

qt 

The first term \jt in rhs of (9.1) generates the Laguerre polynomials. Indeed, it is elementary 
verified that 

, , . 2n + 1 n 2 

On{t) = " , U n [t) = — 
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J — ( 



P n {x- t)P m (x; t)e xt dx + " l ' J mK ' J = h n {t) 5 nm (9.2) 



is a solution of the restricted Toda chain corresponding to the initial condition CQ \t) = 1/1 
These recurrence coefficients correspond to the Laguerre polynomials Lffl(— xt) [13]. 

The second constant term in (9.1) describes adding of a concentrated mass to the measure at 
the endpoint x = of the orthogonality interval. We thus obtain that the polynomials P n (z; t) 
corresponding to (9.1) coincide with the so-called Krall-Laguerre polynomials (see, e.g. [14] for 
details). The measure for the latter is obtained by the adding of (an arbitrary) concentrated 
mass to the point x = of the orthogonality interval for the Laguerre polynomials Lffl(x). In 
our case the orthogonality relation looks as 

<-° P n (0;t)P m (0;t) 
-oo q 
The recurrence coefficients are 

, . . n + 1 2n + 1 

bJt) = h 7 : 

w nt + q (n + l)t + q t 
n 2 (( n -l)t + q)((n+l)t + q) 

Un(t) = ¥ {^TW (9 - 3) 

Remarkably enough that the Krall-Laguerre polynomials belong to a class of 3 families of 
"non-classical" orthogonal polynomials satisfying the ordinary eigenvalue problem for the linear 
differential operator of the 4-th order [14]. 

10. Continued fractions and the Hankel determinants 
connected with the Jacobi elliptic functions 

As we saw, any solution u n (t),b n (t) of the restricted Toda chain is generated by the only 
function c (t). 

The Stieltjes function F(z; t) for the corresponding orthogonal polynomials P n (z; t) is given 
by the Laplace transform 

1 f°° 

F(z;t) = — / c (y + t)e- yz dy (10.1) 
c (t) Jo 

(The factor l/co(t) in front of the integral in (10.1) is needed to provide "conventional" asymp- 
totic behavior F(z;t) = z^ 1 + 0(z~ 2 )). On the other hand, it is well known [8] that any 
Stieltjes function with such asymptotic behavior generates a continued fraction of the Jacobi 
type (so-called J- type continued fraction): 

F(z) = 1 , (10.2) 

z - bo 

Z ~ b l 7 

z - bo - ... 
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where b n ,u n are corresponding recurrence coefficients for the polynomials P n (z;t). Thus we 
can construct explicitly families of continued fractions (10.2) starting from known solution 
u n (t),b n (t) of the Toda chain corresponding to the function co(t) given by (3.4). 

For a special choice of the parameters (3, q, n\ we can obtain families of the Jacobi elliptic 
functions sn(t), cn(t), dn(t) as well as related functions obtained by simplest modular trans- 
forms. Thus we can generalize the Stieltjes results (extended and generalized by Milne [16]) 
who obtained continued fractions corresponding to the Laplace transform of the Jacobi elliptic 
functions in a special case t — of the formula (10.1). 

Moreover, from our results it follows a simple formula for the corresponding Hankel deter- 
minants D n (t) defined by (1.15). 

Indeed, by (1.10) we have 

D n (t) = h (t)h 1 (t)...h n - 1 , n = 1,2,3,... (10.3) 

or, equivalently, 

D n (t) = coit)^- 1 ^- 2 ^) . ..u 2 n _ 2 (t)u n ^(t) (10.4) 
Using explicit formulas (3.2) and (3.1) we find 

h n (t) = e^)+-n!W« , ^ n \\ )w{ \\ P) J q) av (10.5) 

Note that the parameter /x is inessential, because it doesn't contribute to the recurrence co- 
efficients b n (t),u n (t), nevertheless such parameter is convenient when we would like to take 
co(t) coinciding with prescribed Jacobi functions. For the Hankel determinants we have the 
expression 

a(wn(t + 8) + q) .. . . . 

= K « h V exp Aix t + (3) + /i )n) (10.6) 

a{q)a n {w(t + p)) 

where 

/t n = l! 2 2! 2 ...(n-l)!V l(n - 1) (10.7) 

Consider now 3 special cases corresponding to the basic Jacobi elliptic functions sn(t), cn(t), dn(t). 
In all these case we can assume that t\ — e 3 = 1. Indeed, the Jacobi elliptic functions 
sn(t, k), cn(t, k),da(t, k) depend on the modulus 

ei - e 3 ' 

Hence we can pass from given parameters to the scaling parameters je^, i = 1, 2, 3 with some 
nonzero constnat 7. Such transformation doesn't change the Jacobi elliptic functions. Hence 
we can always can assume that e± — e 3 = 1. 
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For the function sa(t) we put w — 1, (5 — q — ujz- We then have the recurrence coefficients 
u 2n (t) = An 2 k 2 (sn 2 (t) - sn 2 (2nt)) 

«w.W = (a, + i)'(*^'(t)- an , ((a ^ + lw ) 

and 

b 2n {t) = (2n + f ) {Z((2n + l)i + itf') + ct/ (2i^)} - 2nZ(2nt) - (4n + l)Z(t) 

Wi(*) = ( 2n + 2)Z((2n + 2)f, fc) - (2n + l)(^((2n + l)i + itf', fc) + in/(2K)) - 
{An + 3)Z(t,k) 

For the function cn(t) we put w — 1, (5 — u 3 , q = uj 2 . We then have the recurrence 
coefficients 



,2,.2/ ...,.2,,, ,/ 2 sn 2 (2nt) \ 
dn 2 (2nt) / 



u 2n (t) = An 2 k 2 -cn 2 (t) + k r 



and 



^ = (2n+ i).(^( t) -*»gg±M) 

, /x / v f w// v ,v sn((2n + l)t, A;)dn((2n + A;) 
MO = (2„ + 1) |Z((2„ + l)t *) - J {2 J + ^ k) 

2nZ(2nt + K,k)- (4n + l)Z(t, k) 



b 2 n+i(t) = (2n + 2)Z((2n + 2)t, k) + K) + 
-(4n + 3)Z(t, fc) 

For the function dn(t) we put w — 1, /3 — u 3 , q = We then have the recurrence 
coefficients 

M2 „(t)=4n 2 f-dn 2 (t)-A;' 2 sn ^\ 



cn 2 (2nt) 

U2n+i(t) = (2n + 1) 2 A; 2 I -cn 2 (t) + k 



2,2 / . 2,,, ,j2 sn 2 ((2n + l)t) 



dn 2 ((2n+ 
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and 



b 2n (t) = (2n + l)Z((2n + l)t + K,k)- 2n\z(2nt, k) - sn ( 2nt ) dn ( 2nt ) 1 _ (4n + 1)z(f> fe) 

J CI1 I ZilXL ) J 

/x / x f„// x ,x sn((2n + 2)t)dn((2n + 2)t)l 
W«) = (2„ + 2) |Z((2„ + 2)t, *) - c , l(( U + l 2 l )t) " } " 

(2n + l)Z(2nt + K,k) - {An + 3)Z(t, k) 



11. Concluding remarks 

The function c (t) given by (3.4) is closely related with simplest solutions of the Lame equation. 
Indeed, consider the Lame equation in the form [2] 

^ = {n(n + l)p(u) + l}y (11.1) 

In case if n is a positive integer one can construct explicit solutions of the Lame equation in 
the form [2] 

<t>{u) = (1L2) 
a n (u) 

where the constants A, ai, . . . , a n can be determined from the Lame equation. In the simplest 
nontrivial case n = 1 we have 

<t>{u) = g (" + fl ) exp(-C(a)«) (11.3) 
a(u) 

where the parameter a is related with the spectral parameter / of the Lame equation by the 
transcendental equation p{a) = I. In case if a ^ u>k, k — 1,2,3 we have the second linearly 
independent solution of the Lame equation in the form 

. . a(u — a) . u . . . . 
4>{u) = K > exp(C(a) M ) (11.4) 
a(u) 

We see that our function co(t) (3.4) coincides with the Lame solution (11.3) when \i\ = — C(o)- 
This means that if the Stieltjes function F(z; t) is the Laplace transform of the solution of 
the Lame equation 

F(z;t)= r + « + e-^Adu (11.5) 

Jo I a(u + t) J 

then corresponding orthogonal polynomials P n (z;t) will have the recurrence coefficients given 
by (3.2) and (3.1) where w — 1, (5 — 0, q — a, //i = — ((a). 
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In [10], [11] some explicit continued fractions connected with the Lame solutions (11.3) 
and (11.2) were announced without any proof or even idea of proof. The authors of [10], [11] 
considered Stieltjes functions of kind of (11.5) but the argument of these function was w = p(a) 
instead of z. This leads to explicit continued fractions which do not resemble presented in the 
present paper. It would be desirable to connect results in [10] and [11] with our ones. 

Another possible generalization consists in passing to the associated polynomials. Indeed, 
we considered here only solutions for the restricted Toda chain, i.e. under the condition u = 0. 
Nevertheless, solutions (3.2), (3.1) can be easily extended to the non-restricted case if one 
replace n with n + c in corresponding formulas, where c is an arbitrary constant not depending 

on t. Then we obtain solution of the nonrestricted Toda chain if c ^ 0,±1,±2, Such 

replacement b n — > b n+c , u n — > u n+c is well known and leads to replacing of the orthogonal 
polynomials P n (z;t) with their c-associated polynomials. Valent already considered [23], [24] 
the c-associated polynomials corresponding to the Stieltjes-Carlitz polynomials. He was able 
to find an explicit orthogonality measure in some special cases. In our case (i.e. for t ^ 0) the 
corresponding analysis seems to be much more complicated. 

Nevertheless, there is a simple special case of the associated polynomials for t ^ which 
leads again to solutions of the restricted Toda chain. 

Indeed, take c = 1 in formulas (3.1) and (3.2) and then put q — 0, w — 1, f3 — 0, jii — 0. 
We obtain the recurrence coefficients 

b n (t) = in + 2)C((n + 2)t) - (n + l)C((n + l)t) - (2n + 3)C(t) (11.6) 

and 

u n (t) = (n + If {pit) - P ((n + l)t)) = ^^^1 (H.7) 

It is seen from (11.7) that u (t) = 0, hence we deal again with a solution of the restricted Toda 
chain. From (11.6) we have 

bo(t) = 2C(2t) - 4C(*) = c (t)/c (t) 

whence 

c (t) = p(t), (11.8) 

where we used the identity [26] 

Thus we obtained that the function Co(i) = p(t) generates another solution of the restricted 
Toda chain described by formulas (11.6), (11.7). This solution was already presented by Chud- 
novsky brothers in [12] The corresponding orthogonal polynomials P n (z; t) seems not to possess 
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positivity property for the Hankel determinants D n (t), hence their measure will not be positive 
on the real axis. 

Nevertheless, if t = 2ui/(N + 2) with some positive integer N we have ttjv(£) = 0. This 
means that polynomials P n (z; t) will be orthogonal on a finite set of points x s : 

N-l 

w s P n (x s ;t)P m (x s ;t) = h n (t) 5 nm , n, m = 0, 1, . . . , N - 1 (11.9) 

where x s are roots of the polynomial Pn(x): 

P n {x s )=0, s = 0,l,...,N-l 

Assume that all the roots x s are simple. Then the discrete weight function w s can be presented 
in the form [8] 

W,= - : (11.10) 

P N ^(x s )P' N (x s ) V ; 

Moreover, for the canonical choice of the parameters e,, i.e. < C2 < e\ we have u n > 0, n — 
1,2, . . . , N — 1 and hence all the weights will be positive w s > 0, s — 0, 1, . . . , N — 1. Finding 
explicit expression for x s and w s is an interesting open problem. 
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